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摘 要 : 针对 极端 频谱 环境 中 强 干 扰 对 全 双 工 系统 性 能 的 影响 ， 提 出 了 一 种 适用 于 频率 捷 变 的 全 双 工 数字 自 干 扰 抑 
制 方法 。 通 过 离线 运算 获取 工作 频带 范围 内 的 误差 补偿 模型 ， 并 利用 该 模型 实时 生成 捷 变 频 点 的 自 干 扰 信道 参数 以 
重建 自 干 扰 信 号 并 完成 数字 对 消 。 数 值 仿真 结果 表明 该 方法 在 系统 工作 频带 范围 内 可 将 自 和 干扰 抑制 到 品 底 以 下 ， 并 
在 整个 频带 有 具有 较 高 一 致 性 和 良好 的 收 全 性 。 
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Digital self-interference cancellation method for frequency agility full duplex system 
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Abstract: To overcome the effects of the intense jam for the full duplex (FD) systems in the extreme Spectrum environment, 
this paper proposed a novel Digital Self-Interference Cancellation (DSIC) method that adapts to the frequency agility full 
duplex systems. The error compensation model of the full working band was calculated offline and could generate the 
self-interference (SI) channel coefficients real-timely to reconstruct the SI signal and accomplish the DSIC. The numerical 
simulation results indicate that the level of the residual SI signal after the DSIC is below the noise floor level in the full band, 
which also demonstrate the consistency and convergence of the proposed method at different frequency points. 
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0 ”引言 .8]、 自 适应 滤波 外 、 辅 助 接收 链 路 的 数字 自 和 干扰 抑制 5.19 等 ， 
到 由 于 需要 重新 估计 自 和 干扰 信号 ， 难 以 满足 同时 同 频 全 双 工 系 
随 着 无 线 通 信 的 飞速 发 展 ， 频 率 资 源 的 需求 日 益 增 大 ， 统 在 动态 频率 捷 变 场景 下 对 自 干 扰 抑制 的 有 效 性 和 实时 性 要 
频谱 竞争 日 趋 激烈 ， 而 传统 的 通信 系统 普遍 存在 频谱 利用 率 求 。 
低 的 问题 。 近 年 来 ， 同 时 同 频 全 双 工 技术 逐渐 引起 人 们 的 关 针对 上 述 问题 , 提出 了 一 种 频率 捷 变 (frequency agility ) 
注 并 取得 了 较 大 的 发 展 ， 该 技术 支持 在 同一 频率 上 同时 进行 的 数字 自 干 扰 抑 制 (digital self-interference cancellation, 
发 射 与 接收 的 通信 和 模式， 可 以 将 频率 利用 率 提 高 一 倍 趾 。 该 DSIC) 方法 ,在 捷 变 工作 频带 范围 内 离线 地 建立 自 干 扰 信 道 
优势 使 得 同时 同 频 全 双 工 技术 成 为 了 学 术 界 和 工业 界 研究 的 自 应 模型 ， 获 取 离 散 的 多 个 频 点 信道 响应 参数 ， 从 而 拟 合 得 
热点 。 到 整个 通 带 的 响应 补偿 模型 。 该 方法 利用 离线 工作 完成 对 全 
同时 同 频 全 双 工 的 核心 问题 在 于 自 干 扰 抑 制 。 扰 抑 通 带 信道 特性 的 估计 ， 在 线 工作 时 实时 地 提供 频率 捷 变 所 对 
制 主 要 从 空域 所、 射频 域 64 和 数字 域名 三 个 维度 展开 ， 通 过 应 的 信道 特性 信息 ,在 通 带 范围 的 任意 频 点 实现 自 干扰 抑 带 
三 个 不 同 作用 域 共 同 作用 的 方式 抑制 自 和 干扰 外 。 在 目前 可 见 满足 同时 同 频 全 双 工 系统 在 频率 捷 变 时 对 自 干 扰 抑 制 的 需 
的 文献 中 ， 研 究 人 员 的 关注 点 主要 集中 于 固定 频 点 、 固 定 带 求 5 
宽 、 外 部 干扰 相对 较 低 的 情况 ， 对 于 极端 频谱 条 件 下 自然 空 让 | 
间或 敌对 环境 中 可 能 存在 的 具有 突 发 性 、 随 机 性 的 强 带 内 干 | 系统 模型 
扰 的 场景 而 言 ， 这 些 干扰 会 对 系统 产生 极 大 的 影响 以 至 于 固 系统 模型 如 图 1 所 示 。 图 1 (a) 中 ， 基 带 产 生发 射 信 号 
定 频 点 、 固 定 带 宽 的 同时 同 频 全 双 工 系统 性 能 下 降 力 至 于 无 x[n] 由 发 射 链 路 处 理 后 得 到 射频 信号 xD ,经 环形 器 后 由 天 线 
法 工作 。 而 频率 捷 变 的 同时 同 频 全 双 工 系统 可 作为 一 个 应 大 发 射 ， 接 收 链 路 接收 到 的 信号 y( 包含 来 自 远 端的 期 望 信号 
的 手段 。 当 出 现 上 述 问题 时 ， 系 统 在 一 定 带宽 范围 内 灵活 调 rlD) 及 环形 器 泄漏 的 扰 信号 *D， 经 过 射频 自 干 扰 抑制 模 
整 自 身 工作 频 点 ， 合 


二 


as 


ES 


理 避 开 强 干扰 频段 。 块 后 得 到 信号 y() ， 到 达 数 字 基 带 后 的 信号 定义 为 yln]， 数 

然而 ， 在 该 工作 模式 下 ， 由 于 工作 的 子 带 发 生变 化 从 而 ” 字 自 干扰 抑制 模块 联合 该 信号 与 发 射 链 路 基带 信号 xj] 生成 
导致 带 内 的 信道 特征 随 之 改变 ， 自 干扰 信号 也 发 生变 化 。 晶 个 估计 信号 Yl] 并 与 yn] 对 消 ， 得 到 zl] 信号 。 图 1 (b) 
前 主流 的 数字 自 干 扰 抑 制 手段 例如 自 和 干扰 信号 估计 与 重建 。” 为 捷 变 模型 ， 信 号 在 射频 链 路 中 传输 时 ， 由 于 切换 频 点 致使 
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信道 响应 改变 , 即 对 应 的 信道 参数 ho 发 生变 化 。 本 文 
同时 同 频 全 双 工 系统 工作 之 前 能 够 获 了 
特性 的 先 验 信息 ， 
出 现时 能 及 时 
应 估计 参数 ， 


保证 自 王 扰 抑制 的 有 效 局 


E 和 实时 性 。 


O ho 外 


Mag 


i fh 


(b) 

(a) 同 时 同 频 全 双 工 系统 模型 (b) 
Fig.1 (a) FD system model; 
同时 同 频 全 双 工 系统 
收发 同 天 线 系统 ， 自 干扰 
及 天 线 反射 0 
性 模型 能 充分 表征 自 干 扰 信 息 ， 
系统 对 自 干 扰 抑 制 的 要 求 。 


Freq 


图 1 


走 变 模型 
(b) Frequency agility model 


言 号 的 主要 来 源 包 才 


| 


期 望 在 
& 整 个 通 带 的 信道 响应 
从 而 建立 误差 补偿 模型 ， 当 系统 捷 变 要 求 
也 通过 误差 补偿 模型 提取 当前 新 频 点 的 信道 响 


扰 主 要 来 源 于 射频 信道 ， 对 于 


笑 环形 器 泄漏 以 
nn 
通过 对 应 的 补偿 模型 能 保证 


考虑 一 个 线性 系统 ,采用 记忆 多 项 式 对 自 干 扰 信 号 建 模 ， 


定义 接收 链 路 加 性 高 斯 白 噪声 (AWGN) 为 uln]， 


到 达 基 带 


的 期 望 信号 (received signal of interest, received SoD) 为 rln], 前 


述 接收 链 路 基带 信号 y'[n] 可 表示 为 
y= hn mtr A (1) 
其 中 : 旋 为 信道 响应 参数 ， 库 sCx% 。 数 字 自 干扰 抑制 模块 


其 从 > 中 减 去 得 到 zl 中] 信和 号 


d= Pha mY hn mt rn 


获取 该 参数 的 估计 值 eC” ， 并 重建 自 干 扰 信 号 3lm] ， 将 


(2) 


上 式 右 侧 前 两 项 为 经 过 DSCI 后 的 残余 


(residual self-interference, residual SI) : 


TResidual SI [a]= haln ml] Thtn m] 

如 何 准确 获取 信道 响应 参数 的 估计 值 是 信 
数字 自 干 扰 抑制 算法 中 的 核心 。 
究 有 较 多 文献 , 其 中 估计 算法 主要 包括 最 小 二 


自 干 扰 信 号 


(3) 
号 估计 与 重建 
对 于 固定 频 点 ， 相 关 算 法 研 
乘 (LS ) 准则 、 


等 : 频率 捷 变 的 全 双 工 数字 自 干 扰 抑 制 方法 


频率 捷 变 范围 为 
点 felfi,fw]， 基 
的 信号 定义 为 


在 离线 工作 ， 


捷 变 动态 场景 模型 


Ph 各 离散 频 点 的 估计 算法 来 
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可 选用 固定 频 点 信道 响应 估计 算法 中 任意 一 种 作为 建立 频率 


[fi,fy]， 对 


于 满足 频率 范围 在 其 
带 信号 x 四 经 过 发 射 与 接收 链 路 后 回 至 


使 用 。 不 妨 设 


中 的 射频 频 


y'iln]= 六 Man —m]+rln]+ ln] 


1 基带 


(4) 


即 无 接收 大 


信号 yiln] 与 信号 


分 别 选 取 这 个 频 点 作为 工作 频 点 习 
系列 不 同 频 点 的 信道 


将 N 个 频 点 归 一 化 至 [0.5.05] 之 间 , 取 互 的 第 m 列 利用 


Xn] 利 | 


D) 。 


望 信号 场景 下 
估计 算法 (如 LS) 得 到 一 组 该 频 点 
的 信道 响应 估计 参数 序列 无 
ho =[h?,hO,.., 
取 子 带 间 隔 为 (fy -fi)/N, 即 产 昌 


rln] 


Pi PD 
hh, ha] 


多 项 式 拟 合 
项 式 拟 合 曲 线 : 


有 .OH 


， 获 得 不 同 频 点 的 入 


gn(fi) 


中 :，g%(f) 为 频 点 处 的 第 m 阶 信道 


响应 估计 参数 值 ，f 为 


=0 。 获 得 的 


(5) 


均匀 分 布 的 N 个 频 点 ， 


EE 复 上 述 工 作 从 而 得 到 一 
响应 估计 参数 太 : 
Bh he io, i = 
[Tho ho hd iD) 
he hh RD 
: 1 (6) 
hi? he ge ho 凶 
[BW BW ,RM) hae) 


和 m 阶 信道 响应 估计 参数 的 多 


(7) 


归 一 化 频率 : fel-0.5,0.5] ，P 为 多 项 式 拟 合 阶 数 ， 4 为 多 
项 式 的 各 阶 系数 。 从 而 ， 计 算得 到 频 点 处 的 所 有 M 阶 信道 响 
应 参数 为 
G(fi) =[goCfi), SIC gm(fi), gui(fi)] (8) 
利用 G(f) 可 快速 计算 获得 捷 变 频率 范围 内 的 任意 频 点 
f 的 信道 响应 估计 值 。 该 算法 按 步 骤 描 述 如 下 : 


a) 离 线 过 程 。 初 始 化 系统 ， 设 置 工作 频带 
子 带 个 数 ， 并 产生 基带 发 射 信号 xn] 。 


ee 


取 接 收 链 路 的 基带 信号 y.'[n] 


的 信道 响应 参数 序列 反 。 


， 并 利用 xj 与 3 


、 子 带 间 隔 、 


c) 离 线 过 程 。 


共有 N 个 子 带 ， 每 个 子 带 对 应 M 个 信道 响 
应 估计 参数 ， 重 复 b) 从 而 获取 所 有 子 带 信道 响应 估计 参数 序 


列 应 =[h0, ho, hd,..., 
M 组 ， 对 每 组 参数 分 别 进行 多 项 式 拟 合 获 得 M 条 拟 合 曲 


CC 。 


d) 在 线 过 程 。 


hwy ， 


将 估计 参数 序列 按 参 数位 


选 定 工 


作 频 率 大 ， 利 月 


线 集合 
自 干 扰 信 号 5 


最 小 均 方 递 归 算 法 (recursive LS，RLS )、 
(MMSE) 准则 等 。 而 对 于 频率 捷 变 的 动态 场景 ， 
见 到 相关 文献 和 工作 报道 。 本 文 提 
自 和 干扰 抑制 方法 ， 解 决 频率 捷 变 下 的 自 干 扰 
2 ， 频率 捷 变 数字 自 干 扰 抑制 


首先 考虑 捷 变 通 带 内 某 一 频 点 的 数字 自 了 


抑制 问题 。 


n 


扰 抑 制 估计 。 


最 小 均 方 误 差 
几乎 没有 
出 了 一 种 频率 捷 变 的 数字 


基于 上 述 方案 ， 
算法 ， 利 用 MATLAB 进行 数值 仿真 分 析 ， 


计算 对 应 频带 的 信道 响应 参数 G(f) , 利用 


有 c) 中 获取 的 参数 


分 为 
线 


该 参数 估计 


] 并 完成 数字 对 消 , 获得 对 消 后 的 信号 区 四 。 
3 ”数值 仿真 与 分 析 


选择 LS 1 


则 作为 离散 频 点 


扰 估 计 


主要 参数 如 表 1 


所 示 。 考 虑 
在 仿真 中 产 
的 FIR 带 
应 。 滤 波 


[ 症 忆 


路 中 不 同 频 点 的 
E 一 个 特殊 的 、 带 内 平坦 度 差 \ 通 带 为 1~ 1.5 GHz 
滤波 器 ， E 变 工作 频带 范 硬 
长 度 为 49 阶 ， 滤 波 特性 


音 道 响应 不 同 


用 来 表征 损 


| 


离线 
LS 准则 


F 程 


期 望 信号 rln]=0 ， 
赴 计 得 到 长 度 为 M 的 参数 序列 


以 及 环 


器 特性 ， 


内 的 信道 响 


如 图 2 所 示 。 


= 


回 到 基带 后 的 信 


号 y,'[n] 通过 


各 和 估计 信号 
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对 消 后 得 到 信号 xs[] 。 自 干 扰 抑 制 效 果 ， 残 余 扰 位 于 噪 底 以 下 ;，(b) 为 数字 
表 1 数值 仿真 参数 干扰 抑制 前 后 接收 期 望 信号 解 调 的 星座 图 ， 通 过 星座 图 证 明 
Table 1 Numerical simulation parameters 了 应 用 该 方法 进行 自 干 扰 抑制 的 有 效 性 。 
参数 说 明 RE 
工作 频带 1~1.5 GHz 有 
子 带 间隔 10 MHz 到 | 必 - 户 =1.26GHe 
子 带 个 数 51 a ee 
发 射 信号 16QAM@50 MHz, 0 dBm 0 
接收 信号 QPSK@20 MHz, -87 dBm 有 
噪声 水 平 -110 dBm 
环形 器 隔离 度 25 dB 已 3 
射频 抑制 50 dB UNUF 
10 10 a pe | 
提 | I | 
本 0 " 5 10 1 20 2 30 
E 5 Polynomial Order P 
号 图 4 多 项 式 拟 合 阶 数 P 与 DCR 的 关系 曲线 
三 -20 Fig.4 DCR vs. Polynomial order P 
-30 -60 
0.3 0.4 0.S 0.6 0.7 0.8 0.9 
Normalized Frequency (x rad/sample) 人 小 
图 2 利用 FIR 带 通 滤波 器 表征 信道 响应 2 
Fig.2 Representing the channel response with FIR bandpass filter EN 
定义 对 消 比 (digital cancellation ratio, DCR ) 为 数字 自 干 3-140 i 
扰 的 抑制 水 平 ， 单 位 为 dB， eve | | 
DCR = Powergso psic -Powerws psic (9) : | | 5 | 
其 中: Powersyore psic 与 Powerwiwr psic 分 别 表 示 数 字 自 干扰 抑制 前 “| i 一 BeforeDSIC 
后 的 信号 功率 。 图 3 为 工作 带 内 不 同 频 点 处 M 取 值 与 对 消 比 30 本 -Rssveasal 
的 关系 曲线 , 在 该 仿真 中 选择 使 各 频 点 对 消 比 在 70 dB 左右 。 2 一 一 一 一 ee 


Digital Cancellation Ratio/ dB 
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Fig.3 DCR vs. Estimation order M 


通过 离线 过 程 获得 了 了 


[ 作 


频带 中 均匀 分 布 的 51 


的 信道 响应 参数 序列 ， 
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获取 长 


度 为 M=40 的 拟 合 


线 集合 G(f) 


按照 前 述 规则 进行 多 项 式 拟 合 以 
。 图 4 为 多 项 式 扫 


阶 数 己 与 对 消 比 之 


间 的 关系 


线 ， 选 择 与 


图 3 中 测试 频 点 相 


(a) Baseband 


频率 ， 纵 多 
个 频 点 


全 
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位 置 效 果 很 差 。 
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(9) 射 频 工 作 频 点 f=1.46GHz 处 基带 频谱 
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spectrum for Radio frequency point: f=1.46 GHz 
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(b) Constellation results before and after DSIC 
图 5 仿真 结果 


Fig.5 Simulation results 


工作 全 频带 仿真 结果 如 图 6 所 示 ， 
E 标 分 别 为 EVM 
文 提出 的 方法 可 将 EVM 控制 在 8% 左 右 ， 
70 dB 左右 , 与 离线 过 程 吕 
定 参 数 下 的 DSIC 仅 在 部 分 点 的 抑制 ? 
6 证 明了 本 文 所 提 


多 


P 直接 估计 得 到 


其 中 横 坐 标 为 归 一 化 


与 DCR。 由 图 6 (a) 可 见 利 


本 
自 干 扰 抑制 量 在 
上 效果 相 一 致 ; 
而 在 其 他 


的 抑 表 
效果 较 好 ， 


出 自 


差 10 M 的 新 的 测试 频 点 进行 测试 ， 从 图 中 可 以 看 出 ， 当 拟 范围 内 的 一 致 性 与 有 效 性 。 
合 阶 数 为 P=21 阶 时 , 新 的 测试 频 点 均 能 达到 70 dB 左右 的 自 十 上 所 :五 

8 4 结束 语 
干扰 抑制 效果 。 

在 线 运行 时 ， 在 信号 中 加 入 接收 的 期 望 信号 (Received 本 文 针 对 频率 捷 变 场 景 下 信道 响应 改变 致使 数字 自 干 扰 
SoI)， 利 用 离线 过 程 生 成 的 多 项 式 提取 估计 参数 用 于 数字 自 抑制 过 程 中 重建 自 干 扰 信 号 误差 较 大 而 引起 的 自 干 扰 抑 制 性 
干扰 抑制 ， 仿 真 结果 如 图 5 所 示 ， 图 中 ，(a) 为 射频 工作 频 能 降低 的 问题 ， 提 出 了 一 种 离线 建立 工作 全 频带 内 的 响应 补 
点 f=1.46GHz 处 的 基带 频谱 , 通过 该 方法 获取 了 70dB 左右 的 “，” 偿 模 型 ,并 在 频率 捷 变 过 程 中 利用 该 模型 实时 提取 补偿 参数 ， 
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频率 捷 变 时 参数 的 
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建 自 和 干扰 信号 的 方法 对 自 干 扰 进 行 抑制 的 方法 。 该 方法 在 
离线 过 程 完 成 了 响应 模型 的 建立 从 而 避免 了 在 线 运行 过 程 中 
生计 ， 满 足 动态 场景 下 自 干 扰 抑制 对 
时 效 性 的 要 求 。 数 值 仿真 结果 显示 该 方法 能 实现 通 带 范围 内 


的 进一步 优化 ， 并 进行 实际 系统 的 实验 验证 等 。 
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(b) 工 作 全 频带 各 频 点 DCR 曲线 
(b) DCR of full working band 
图 6 工作 全 频带 仿真 结 
Fig.6 Working full band simulation results 
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